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CONTRIBUTION TO METHODS FOR CALCULATING THE FLOW 


ABOUT THIN LIFTING WINGS AT TRANSONIC SPEEDS - 
ANALYTIC EXPRESSIONS FOR THE FAR FIELD 

By E. B. Klunker 
Langley Research Center 

SUMMARY 

The problem of determining the small-disturbance flow about two-dimensional air- 
foils at transonic speeds has been successfully treated by Murman and Cole by the 
process of matching a numerical solution of the near field to analytic expressions for the 
far field. The three-dimensional problem, it would appear, can be treated in a similar 
way with the aid of algorithms adapted to high-speed and high-capacity computers. The 
far-field potential for both lifting and nonlifting three-dimensional wings at transonic 
speeds is developed herein for a subsonic free stream. This potential could be used for 
a three-dimensional -wing computation similar to the computation made by Murman and 
Cole for the two-dimensional wing. 


INTRODUCTION 

The problem of calculating the flow about two-dimensional airfoils at transonic 
speeds with a subsonic free -stream Mach number has been successfully treated by 
several authors using both time-dependent and steady-flow calculation techniques. The 
elliptic nature of the problem, which requires that the disturbances vanish at infinity, 
leads to the use of a large computational network together with large demands on computer 
storage unless special provision is made for satisfying the boundary condition at infinity. 
The numerical solution of a three-dimensional problem is even more demanding on 
computer storage and time requirements, and it is essential to restrict the computational 
network insofar as possible in order to obtain solutions with present-generation computers 
in a reasonable time. 

Several techniques have been developed for satisfying the far-field boundary condi- 
tions in the finite-difference computations of the flow about two-dimensional airfoils. 

Sills (ref. 1) has taken the simplest approach by simply extending the computational net- 
work sufficiently far from the airfoil to satisfy effectively the boundary conditions in the 
far field. A large computational network is required with this approach, particularly for 



lifting airfoils at transonic speeds, since the disturbances created by the body decay 
slowly in the lateral direction. Magnus and Gallaher (ref. 2) limited the computational 
grid by transforming the far-field region into a finite region by an inversion. Matching 
the computational mesh of the near field with that of the far field presents a difficulty with 
this method. The transformation of the airfoil into a circle and a subsequent inversion 
has been used by Sells (ref. 3) and others to restrict the computational network to the 
interior of a circle. This technique leads to a nearly ideal computational system for two- 
dimensional flows, but the extension to three-dimensional flows is not evident. The 
method developed by Murman and Cole (ref. 4) provides a means for limiting the compu- 
tational network and satisfying the far-field boundary conditions. Their approach is to 
match the numerical solution of the near field to an analytic representation of the far field 
to satisfy, in effect, the far-field boundary conditions. This technique can be extended 
to three-dimensional wings, for which an economy of grid points is a prime concern. 

The Murman and Cole technique requires an analytic expression for the asymptotic 
form of the potential at infinity. A rigorous development of the asymptotic form of sub- 
subsonic two-dimensional flow solutions has been given by Ludford (ref. 5) as a series. 
The leading terms for the far-field behavior for the two-dimensional airfoil at transonic 
speeds are given by an integral representation for the nonlifting case in reference 4 
and for the lifting case in reference 6. The purpose of the present paper is to develop 
the leading terms of the far-field representation for both lifting and nonlifting three- 
dimensional wings appropriate to a transonic flow with a subsonic free stream. The 
development is in the form of an integral representation similar to that of reference 4. 

The two-dimensional results are presented as a special case. 

SYMBOLS 


b wing span 

wing thickness distribution 
i,j,k unit vectors along X-, Y-, and Z-axes 

L Laplacian operator 

M w free-stream Mach number 

n inward normal direction 

O order symbol 


2 



R 


magnitude of difference vector defined by equation (6) or (14) 


S 

u,v,w 

V 

X,Y,Z 
x,y,z 
x,r ,n 

0 = fl - M~ a 
r 

v 

A 

<#> 

□ 

Subscripts: 

D 

1 


surface 

disturbance velocity component in x-, y-, and z-direction, respectively 

J 

volume 

Cartesian axes 
Cartesian coordinates 

cylindrical coordinates 
2 

circulation 

ratio of specific heats 

denotes difference between value on upper and lower surface of wing 
dummy variables 
disturbance velocity potential 
solution to Laplace equation 

angles between shock normal and coordinate axes 
denotes jump values across shock 

shock discontinuity 
interior region 


P 


singular point (x’,y',z') 



V 


trailing vortex sheet 


W wing 

1,2 upstream and downstream side of shock, respectively 

A prime indicates a transformed space in which the coordinates transverse to the 
stream direction are stretched by the factor j3. Single and double subscripts on 0 
and \p denote first and second derivatives with respect to the indicated variables. 

ANALYSIS 

In the following development the transonic small-disturbance differential equation 
is put in the form of an integral equation (with the use of Green's theorem) which is 
suitable for representing the flow about three-dimensional wings at transonic speeds with 
a subsonic free stream. The leading terms of the asymptotic form of the equation at 
infinity are taken as the far -field representation. 

Transonic Differential Equation 

Let x, y, and z be Cartesian coordinates and let <p be the disturbance velocity 
potential that is the perturbation about the uniform free stream. The free stream is 
taken in the x-direction and the wing is in the XY-plane. The small-disturbance potential- 
flow equation applicable to transonic flows is 

+ «Vy + ^zz = i M (1) 

where = 1 - The disturbance velocity components are u = <p x , v = $y, 

and w = <£ z . The left-hand side of equation (1) becomes the Laplacian operator 
for p? > q W ith the transformation x' = x, y' = /3y, and z' = /3z. In terms of x', y', 
and z', 

L(0) = 4>x'x' + 0y'y' + ^z'z' = ^2 Sx 7 (^ x ' ) ^ 

In the next section an integral equation for the velocity potential 4> is developed from 
equation (2) and subsequently transformed back to the unprimed coordinates. 

Integral Equation 

The symmetric form of Green's theorem for the two functions <£ and 4/ which 
are twice continuously differentiable in the volume V bounded by the surface S is 


4 


(3) 


j^LW.) - »LW) dV » dS' 

where n is the inward normal direction to S and where L is the Laplacian operator. 
Let 


L(^) =^ x , x , + ^ y ,y,+ ^ Z . Z .= 0 

The fundamental solution of equation (4) is 



where 



n 2 + (y 


V') 2 + (z' 



(4) 

(5) 

( 6 ) 


This choice of 4* leads to source and doublet distributions in the integrals over the wing 
surface. The boundary surface S (fig. 1) includes (1) the surface at infinity, (2) the 



Figure 1.- Geometric surface comprising boundary surface S. 
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surface Sp surrounding the singular point (x',y',z'), (3) the surface Sp around any 
shock discontinuity, (4) the wing surface and (5) the trailing vortex sheet Sy. The 

integral over the surface at infinity, with ifs given by equation (5), vanishes under the 
condition 4> ~ R“ , for e > 0, everywhere at infinity except in a region surrounding the 
trailing vortex sheet where <p is required to be bounded and antisymmetric only in z. 
The far -field expression for <p will be shown to have these properties and hence the 
integration over the surface at infinity vanishes and need not be considered further. 

In the following development the potential at the field point (x',y',z’), which is 
obtained by integration over the surface Sp, is related to the remaining surface and 
volume integrals with the use of equation (3). The integrals over any shock discontinuity 
combine, so the integrand vanishes and thus the combined integral makes no contribution 
to <p. The integration over the trailing vortex sheet is transformed through an integra- 
tion by parts so that the potential at the field point is given in terms of only surface inte- 
grals over the wing, which represent source and vortex distributions, and a volume inte- 
gration throughout the field. 

Integral over volume V .- Since L(i^) = o, the only contribution to the left-hand 
side of equation (3) comes from the term \]/L,((p). The reduction is made with the use of 
the right-hand side of equation (2) for L($) together with an integration by parts in the 
x' -direction. The integration by parts leads to integrals over the bounding surface of the 
volume V, all of which vanish except those over any shock discontinuity. Thus, 


1 


i//L(0)dV' = 


y + 1 
2/3 2 



(7) 


where the jump across a shock surface is 



" U 1 


2 


with the subscripts 1 and 2 referring to the upstream and downstream values, respec- 
tively, at the shock-discontinuity surface Sp. The upstream unit vector normal to the 
shock has been taken as 


n = i cos + j cos ^ 2 ' + k cos ^ 3 ' 


and in equation (7) 

cos Si’dS' = sgn^cos Q^d?]’ d = -dry' d£’ 
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Integral over surface surrounding singular point Sp.- The field point (x',y',z') is a 
singularity of the equations for R = 0 and must be excluded from the region of integra- 
tion. The surface surrounding the field point is taken as the sphere R = R 0 . Then 


8 £_ _ 1 
8n 4 ttR 0 2 

and the element of area is dS' = R 0 2 da>, where u> is the solid angle. Integration over 
the surface Sp gives 



ip I^AdS' = lim ^ r 4 ", 

9 «y R n _o 4,r J o I 


1_ 


R 2 



Ro **0 j 


Ro 2 d 


CO 


= -cp(x\ y',z') 


(8) 


Integral over surface of discontinuity Sp.- The quantities xp, J^-, and (p are 
continuous across any shock surface and the normal velocity is discontinuous there. 

Thus, these surfaces must be taken as boundaries of the region of integration in order to 
satisfy the requirements of Green's theorem The integral over Sjj becomes 



(9) 


where 


denotes the jump in 


across the shock and 


upstream normal to Sq. Equation (9) will be used with 

[ff] = -[VJcos V - [V] COS " p?'] COS 

Integral over wing surface S 


> w .- For a thin wing — = ± ^ 


0 ^- is in the direction of the 


expressed in the form 

( 10 ) 

on the wing surface. 


9 


It is convenient to employ the notation A</> and A to designate the value on the 
upper surface minus the value on the lower surface. Thus, A (p is the jump in potential 
across the wing surface and A is related to the wing thickness distribution. The 

d£' 

integral over the wing surface then becomes 


- * !£) 1 as ' = ^ - +* !$) ds ' (I1 » 

Integral over trailing vortex sheet Sy.- For lifting wings there is a trailing vortex 
sheet Sy which extends downstream from the wing trailing edge and lies in the 
X’Y' -plane in the small -disturbance approximation. The quantities ip, and 
are continuous on Sy whereas <p is discontinuous there. Thus, 
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I 



H A + ds ’ 


( 12 ) 




Integral equation for velocity potential .- Equations (7) to (12) provide the relations 
required to relate the potential at the field point to the surface integrals over the wing, 
the trailing vortex sheet, and any shock discontinuities as well as to a volume integral 
throughout the space. Substitution into equation (3), together with a transformation to 
the unprimed variables, gives 


* = A * ■ * iw ) dS * ^ f ** dS * ^ ^ 

[u2]} c ° s + [jrjcos ^2 + [VJcos 9 


3 ) dS (13) 


where 


^ = 


4irR 

with R now defined as 


R 


•{- 


*) 2 + 


P 2 (y - v ) 2 


/3 2 (z 



(14) 


and where the velocity components are u = 0^, v = <p and w = <p^. The relations 
cos dS’ = p^cos 8 ^ dS 
cos $ 2 ' dS' = p cos ^2 dS 
cos £ 3 ' dS' = p cos 83 dS 

have been used in the transformation of the surface integrals from the primed to unprimed 
coordinates. 

The direct expansion of the shock polar equation for weak jumps gives the equation 
js 2 DO - 14 a [ u 2]j M + 002 + 002 = „ (15) 


8 



which may be regarded as a small-disturbance approximation to the shock jump condi- 
tions. Applying the divergence theorem to the irrotationality conditions gives 


JjQcos - QTJcos ^2 = 0 
[vf]cos - QT] cos ^3 = 0 


cos ^3 = 0 

and these relations combined with equation (15) give 


races s 2 - [x i 


2 M - ^[«9}cob % 


+ Cvjcos ^2 + ra c °s ^3 = 0 


(16) 


as a second form of the shock jump conditions for weak jumps. The left-hand side of 
equation (16) is a factor in the integrand of the integral over Sj) in equation (13) and, 
consequently, that integral does not contribute to the potential <p. The two forms of the 
shock jump conditions are equivalent to those given in reference 4. 


The terms in equation (13) with integrands A <f> correspond to doublets, and the 

region of integration is over both the wing and the trailing vortex sheet. These integrals 
can be transformed to one extending only over the wing surface through integration by 
parts in the ij -direction and by making use of the fact that A <j> is a function only of r\ 
on Sy. The potential becomes 


4 > = 



Au 

(y - fj ) 2 + z 2 



Aw dS + 



u fyg dV (17) 


The first term on the right-hand side of equation (17) corresponds to a vortex distribution 
and represents the lifting effects; the second term corresponds to a source distribution 
whose strength is related to the wing thickness distribution; and the third term, which 
arises from the nonlinearity of the flow equations, has the form of a doublet with its axis 
in the stream direction and with a strength given by the local value of u 2 and influences 
the results for both lifting and nonlifting wings. 

The corresponding relation for two-dimensional flow is found from equation (17) by 
regarding Au, u, and Aw to be functions only of the argument £, and integrating 
on 77 from -00 to +°o. Thus, for a two-dimensional flow, 
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( 18 ) 


<t> = ± f Au ft 

2ir Jq |^2 


sgn(z) + tan 


-1 


+ (3*z 


2„2 


>W 


J 


4tt0 


H" 2 


S ^ i )} d{ + 253^ awlog { <x -* )2 

(x - £)dS 

(x - £) 2 + /3 2 (z - S) 2 


where C designates the airfoil section and S is the two-dimensional region exterior 
to the airfoil. The term sgn(z) in the first integral is required to account for the jump 
in potential downstream of the trailing edge on z = 0 with the principal value taken for 
the arc tangent. Equations (17) and (18) provide the basis for obtaining the asymptotic 
form of the potential in the far field for the three-dimensional and two-dimensional wings, 
respectively. 


Far-Field Representation 

In the far field the distance from the wing is large compared with the dimensions 
of the wing. Each integral in the expression for the potential given by equation (17) or (18) 
can be approximated and simplified at large distances from the wing on this basis. 


Lift integral.- There are two regions of the far field to consider for the approxima- 


tion of the integral associated with the lift: (1) the region where the distance from the 

wing (x2 + y2 

' ' ' y 


5 2 ) 1/2 i ar g e w tth both y and z exterior to a region 0(b) 

surrounding the trailing vortex sheet and (2) the region where only x is large and 
and z lie in a neighborhood of the trailing vortex sheet. 

For the first region, where (x2 + y 2 + z 2)*/ 2 is large compared with the wing 
dimension, 


x - £ — x 

y - 17- y 



and the contribution from the integral associated with the lift in equation (17) becomes 


hift 


1 z 
4 tt y 2 + z 2 



( 19 ) 


For the second region (in the far field in the neighborhood of the trailing vortex 
sheet), x »b and both y and z are 0(b). The term in parentheses in the first 




integral in equation (17) approaches a value of 2 as x 
integral associated with the lift is 


^lift 


\ 

2ir J_ 


pb/2 

-b/2 (y - 7?) 2 + z 2 


IM. 


d rj 


oo, and the contribution of the 


(20) 


where the circulation r is 


r « * 1 


Au(£,? 7 )d£ 


and the integration is over the chord. On the trailing vortex sheet, the only contribution 
from equation (20) comes when the denominator approaches zero. This limiting form is 
readily evaluated as 


^lift 

*lift 


~ T(y) 
2 

2 



( 21 ) 


The two approximations given by equations (19) and (20) are in agreement outside a region 
where y is 0(b) and z is 0(b). 


Thickness integral .- In equation (17) the integral with the integrand Aw is 
related to the wing thickness distribution through the surface boundary condition 


Aw = 2F^ (£ ,rj) 

where F (£, 77 ) is the wing thickness distribution. Integrating by parts and letting £ 
and ?? approach zero gives the leading term of the far-field representation for the non- 
lifting part of the potential as 

^thickness ~ "^3 F(M)dS (22) 

which is of higher order than equation (19). 

Volume integral .- The contribution to the far-field potential from the volume inte- 
gral is shown to be of higher order than that associated with the lift. The volume integral 
in equation (17) extends over an infinite region; however, this integral can be limited to a 
finite interior region Vj with little error since u decays rapidly far from the wing. 

Differentiation of equation (17) shows that u ~ -ts- for large R. Then the contribution to 

R^ 

the potential from the region exterior to Vj is 


11 
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and the contribution from the neighborhood of the field point is of even higher order. 


The order of magnitude of the far -field potential due to the integral over the 
interior region Vj can be estimated with the use of the Schwarz inequality 




^ 2 dV 


p 

The integral of can be obtained in closed form in cylindrical coordinates and the 

result can then be expanded in powers of ^ where r is related to the region of inte- 

rt 

gration as shown in figure 2. The region of integration is the annular sector defined by 


r-rS=r = r + ? 



-f^^r 


(y.z) 



Figure 2.- Polar coordinates in YZ-plane. 




where r = fify* + z 2 and /j. is the polar angle (fig. 2). The Schwarz inequality gives 




( 23 ) 


where 


* V, ^ 


dV 


and the volume of the interior region is 
Vj = ( 2 r ) 3 

Equation (23) yields 


j* u 2 i^ dV ~ O 


Thus the dominant contribution to the far-field potential from the volume integral is 
<p ~ ^ ■ * * ^ u 2 i^ dV 

which is Of -^ 5 - 
\R 3 

Veloc ity potential for the far field .- The leading terms of the potential for the far 
field are 


V Mds+ ^V Vv 

where the various forms of the contribution due to lift 0^ are given by equations (19) 
to (21). The terms associated with the lift largely dominate the far-field potential. 

These terms provide the proper jump across the trailing vortex sheet and they decay like 
the reciprocal of the distance in the regions above and below the wing as well as far down- 
stream, whereas the contributions from the thickness and the volume integrals are 

It should be noted that 4>, as required by Green's theorem, has the proper form 

far from the wing (<f> ~R~ e , for e > 0 , and <p is bounded and antisymmetric in z near 
the trailing vortex sheet) to ensure that the integral over the surface at infinity vanishes. 
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The far-field representation for the two-dimensional airfoil is found from equa- 
tion (18) to be 


0 


2 v 



sgn(z) + tan 


-1 



+ -1 * 

X 2 + j 32 z 2 


\ nm 

’c 


+ 


y + 1 

4ir/3 


? u 2 (x_- g) d g . . . 

’Sj (x - S) 2 + <3 2 (z - 0 2 


where C designates the airfoil section and Sj is the two-dimensional equivalent of 
region Vj. The two-dimensional potential is singular for = 1; the free-stream 
Mach number can approach 1 only as the thickness ratio t simultaneously approaches 0 
so that the transonic similarity parameter /3 2 /t 2 /3 remains of order 1. 

Application o f the far-field equation.- The far-field representation developed herein 
would be used for a three-dimensional -wing computation in the same manner as the 
corresponding two-dimensional far-field representation was used in the numerical calcu- 
lations of Murman and Cole (ref. 4) and Krupp and Murman (ref. 7). A numerical solu- 
tion of equation (1) would be made within a limited region surrounding the wing and the 
values of the potential would be matched to the far-field values on the outer boundary of 
the computational grid to satisfy, in effect, the boundary conditions at infinity. Only that 
part of the potential which is associated with the lift would be used for the far-field 
representation for lifting wings since the other terms are of higher order and may be 
neglected. The far-field value of the potential associated with the lift as given by equa- 
tion (19) would be used everywhere except in some neighborhood of the trailing vortex 
sheet far downstream of the wing where equation (20) would be used. The region where 
the different forms of the equation for (see eqs. (19) to (21)) merge would have to 

be established by trial. The integrands in the far-field representation contain the 
unknown velocity component u and consequently have to be evaluated during the finite- 
difference computation. Thus an iterative procedure is required in the numerical 
process. 


CONCLUDING REMARKS 

A far-field representation for the transonic flow about lifting wings is developed 
for the case where the free-stream Mach number is subsonic. The term associated with 
the lift is shown to dominate the far-field solution. This representation provides a means 
for both limiting the computational network for a finite -difference computation (thereby 
limiting the computer storage requirements) and satisfying the requirement of vanishing 


I 
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disturbances at infinity. A finite -difference computation would be made within a limited 
region surrounding the wing and the values of the potential would be matched on the outer 
boundary of the computational grid to satisfy, in effect, the boundary conditions at infinity. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 24, 1971. 


REFERENCES 

1. Sills, J. A.: Computation of Inviscid, Two-Dimensional, Transonic Flows Using a 

Time-Dependent Finite-Difference Method. ERR-FW-806, Gen. Dyn./Convair, 
Dec. 18, 1968. 

2. Magnus, Richard J.; and Gallaher, William H.: Flow Over Airfoils in the Transonic 

Regime — Computer Programs AFFDL-TR-70-16, Vol. II, U.S. Air Force, 

Mar. 1, 1970. 

3. Sells, C. C. L: Plane Subcritical Flow Past a Lifting Aerofoil. Proc. Roy. Soc. 

(London), ser. A, vol. 308, no. 1494, Jan. 14, 1969, pp. 377-401. 

4. Murman, Earll M.; and Cole, Julian D.: Calculation of Plane Steady Transonic Flows. 

AIAA J., vol. 9, no. 1, Jan. 1971, pp. 114-121. 

5. Ludford, G. S. S.: The Behavior at Infinity of the Potential Function of a Two Dimen- 

sional Subsonic Compressible Flow. J. Math. Phys., vol. XXX, no. 3, Oct. 1951, 
pp. 117-130. 

6. Cole, Julian D.: Twenty Years of Transonic Flow. Dl-82-0878, Flight Sci. Lab., 

Boeing Sci. Res. Lab., July 1969. 

7. Krupp, J. A.; and Murman, E. M.: The Numerical Calculation of Steady Transonic 

Flows Past Thin Lifting Airfoils and Slender Bodies. AIAA Paper No. 71-566, 

June 1971. 


NASA-Langley, 1971 


12 L-7917 


15 



NATIONAL AERONAUTICS AND SPACE ADM I STRATI ON 
WASHINGTON, D.C. 20546 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


off.cal busin ess FIRST CLASS MAIL 

PENALTY FOR PRIVATE USE S300 



018 001 Cl U 12 711112 S00903DS 
DEPT OF THE AIR FORCE 
AF WEAPONS LAB <AFSC) 

TFCH L IBRARY/WLOL/ 

ATTN: F LOU BOWMAN, CHIEF 
KIRTLAND AFB NM 87117 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual ) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the at?)iosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 

NASA S CIE NTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech BriefSj 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




